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Abstract 

We analyze several problems of Optimal Transport Theory in the 
setting of Ergodic Theory. In a certain class of problems we consider 
questions in Ergodic Transport which are generalizations of the ones 
in Ergodic Optimization. 

Another class of problems is the following: suppose a is the shift 
acting on Bernoulli space X = {0, 1} N , and, consider a fixed contin- 
uous cost function c : X x X — > M. Denote by IT the set of all Borel 
probabilities tt on X x X, such that, both its x and y marginal are a- 
invariant probabilities. We are interested in the optimal plan tt which 
minimizes f c dir among the probabilities on IT 

We show, among other things, the analogous Kantorovich Dual- 
ity Theorem. We also analyze uniqueness of the optimal plan under 
generic assumptions on c. We investigate the existence of a dual pair of 
Lipschitz functions which realizes the present dual Kantorovich prob- 
lem under the assumption that the cost is Lipschitz continuous. For 
continuous costs c the corresponding results in the Classical Transport 
Theory and in Ergodic Transport Theory can be, eventually, different. 

We also consider the problem of approximating the optimal plan 
tt by convex combinations of plans such that the support projects in 
periodic orbits. 



1 Introduction 



For a compact metric space X, we denote P{X) the set of probabilities 
acting on the Borel sigma-algebra B(X). C(X) denotes the set of continuous 
functions on X taking real values. 
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We denote by a the shift acting on {1,2,.., d} N and by a the shift acting 
on {1, 2, .., d} z . Some of our results apply to more general cases where one can 
consider a continuous transformation defined on any compact metric space X. 
Anyway, the reader can take {1, 2, .., d} z = XxY = {1, 2, .., d} N x {1, 2, .., d} N 
as our favorite toy model. 

We consider a continuous cost function c : X x Y — > IR, where X and Y 
are compact metric spaces. 

The Classical Transport Problem consider probabilities tt on P(X x Y) 
and the minimization of J c(x,y)d7r(x,y) under the hypothesis that the y- 
marginal of tt is a fixed probability v and x-marginal of tt is a fixed probability 
fi. A probability tt which minimizes such integral is called an optimal plan 

We want to analyze a different class of problem where in some way the 
restriction to invariant probabilities [22] appears in some form. We present 
several different settings. 

As a motivation one can ask: given the 2-Wasserstein metric W on the 
space of probabilities on X = {1, 2, .., <i} N , and a certain fixed probability fi, 
which is not invariant for the shift o : {l,2,..,<i} N — > {1, 2, .., d} n , charac- 
terize the closest a-invariant probability v to fi. In other words, we can be 
interested in finding a a-invariant probability v which minimizes the value 
W(fi, v) for a fixed fi. In this case we are taking X = Y. What can be said 
about optimal transport plans, duality, etc? 

As a generalization of this problem one can consider a continuous cost 
c(x, y), where c : {1, 2, .., d} N x {1, 2, .., <i} N — > M, and ask about the properties 
of the plan tt on {1, 2, .., d} N x {1, 2, .., <i} N which minimize J c(x, y)dir(x, y) 
under the hypothesis that the y-marginal of tt is a variable cr-invariant prob- 
ability and the x-marginal of tt is a fixed probability /i. 

We note that a plain with this marginals properties is characterized by: 

f / f(x) dir(x, y) = J f(x) dn(x) for any / e C(X) ,^ 
\ J g(y)dTr(x,y) = J g(a(y))dTr(x,y) for anyg e CiY) 

We will show in section [2] the following: 

Theorem 1 (Kantorovich duality). Consider a compact metric space X 
and Y = {1,2,.., d} N . Consider a fixed fi G P(X), and a fixed continuous cost 
function c : X x Y — > M + . Define H(fi, er) as the set of all Borel probabilities 
tt G P(X x Y) satisfying (QJ]. Define <3> c as the set of all pair of continuous 
functions G C(X) x C(Y) which satisfy: 

<f(x)+i[)(y)-(i/;oa)(y)<c(x,y), W(x,y)eXxY (2) 

Then, 
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I) 



inf 

n(/i,<j) 



/ 



cdn = sup 

(t£,V>)e$ c 



/ 



(3) 



Moreover, the infimum in the left hand side is attained. 

II) If c is a Lipschitz continuous function, then, there exist Lipschitz (p 
and ip which are admissible realizers of the supremum. 

Any pair (p and ip satisfying (T5]) is called admissible. Any it realizing 
the infimum in ([3]) will be called an optimal plan, and its y-projection an 
optimal invariant probability solution for c and fi. Moreover, p and ip are 
called an optimal dual pair if they realize the maximum of the right hand 
side expression. It is possible that does not exists an optimal dual pair (see 
remark bellow). 

The following criteria is quite useful. 

Slackness condition [27] [22]: suppose for all (x,y) in the support of 
vr G n(yU, a) we have that 



for some admissible p and xp satisfying p + ip — {ip o a) < c, then tc is an 
optimal plan and p and ip is an optimal dual pair. 

In recent years several results in the so called Ergodic Optimization The- 
ory were obtained [E][6][I3[2][23][3][l0]p. We will show that the above 
kind of Ergodic Transport problem contains as a particular case this other 
theory. The subaction which possess properties of minimality described in 
[5] and [TT] can be seen as a version of Kantorovich duality. 

In the below remark we show that there are conceptual differences in 
the kind of analogous results we can get in the Classical and in the Ergodic 
Transport Theory. 

Remark 2. We point out that if fi is a Dirac delta in a point xq, then the 
cost c(xo,y) just depends on y. In this way if we denote A(y) = c(xo,y) 
we get that the above problem is the classical one in Ergodic Optimization, 
where one is interested in minimizing J Adv among invariant probabilities 
v. There is no big difference in this theory if one consider maximization 
instead of minimization. The function ip above corresponds to the concept of 
subaction and the number p(xo) is equal to min{J A(y) dv(y) : v is invariant} 
W$' W$' W$' llM- ft i s known that for a C°-generic continuous potential A 
does not exist a continuous subaction %4jj. For the Classical Transport problem 
in compact spaces there exists continuous realizers for the dual problem when 
c is continuous W8$ . This shows that there are non trivial differences (at 
least for a C° cost function c) between the Classical and the Ergodic transport 



<p(x) + ip(y) - (ip o a){y) = c(x, y), 
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setting. It is known the existence of a calibrated Holder subaction for a Holder 
potential A. The item II on the above theorem is the correspondent result on 
the present setting. The expression ip(x) + i>(o~(y)) — ip(y) < c(x,y) can be 
for some people more natural. This can be also obtained by replacing if) by 

The next example shows that the Ergodic Transport problem can not be 
derived in an easy way from Ergodic Optimization properties. 

We denote by (ai...a n )°° the periodic point {ax...a n ai...a n a\...) in {0, 1} N . 

Example 3. Consider X = {xo,xi}, \i = h(8 X0 + 6 X1 ), Y = {0, 1} N , and a 
cost function c defined on X x Y , satisfying the following proprieties: 

1) c(ar , (01)°°) = 0, c(ar , (10)°°) = 1, c(x Q ,W°) = 1/4, c(x ,y) > 0, if 
y?(pi)°°. 

2) 0(^,(01)°°) = 1, 0(^,(10)°°) = 0, 0(^,1°°) = 1/4, c( Xl ,y) > 0, if 
y^{10)°°. 

Assume c is Lipschitz continuous. 
Note that, as an example, we can take 

c(x , y) = d\y, (01)°°), c(x u y) = d\y, (10)°°). 

We observe that the measure v = o(^(oi)°° + ^(io)°°) ^ s n °t a minimizing 
measure for either of the potentials A (y) := c(x ,y), or, A 1 (y) := c{x\,y). 
By the other hand, the unique optimal plan is given by 

71 = ^(fyso.COl) 00 ) + ^1,(10)°°)), 

which projects on fi and v. 

We will also show in section [3] that generically on c the optimal plan is 
unique. 

In another kind of problem one can ask: given a continuous cost c(x,y), 
c : {1, 2, .., d} n x {1, 2, .., d} N — > R, what are the properties of the probability 
7r on {l,2,..,<i} N x {1,2, ..,d} N which minimize J c(x,y)dn(x,y) under the 
hypothesis that the y-marginal of tt is a variable invariant probability v 
and the x-marginal of it is a variable invariant probability /i? Under what 
assumptions on c we get that the optimal plan n is invariant for a? 

We will present now formal definitions of the second class of problems. 

Here we fix compact metric spaces X and Y and continuous transforma- 
tions 

Ti : X x Y -)• X, T 2 :XxY^Y, 
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such that, T : XxY XxY, given by T = (Ti, T 2 ), defines a transformation 
of X x F to itself. Let IT(T) the set of Borel probability measures tt in X x F, 
such that, for any / : X ->■ R, # : F ->■ R: 

f(x)dn(x,y) = J f{T 1 {x,y))d'K{x,y), 

and 

g(y)dn(x,y) = / g(T 2 (x,y)) dn(x,y). 



The set of such 7r is called the set of admissible plans. 
Note that any T-invariant measure in X x F (which exists because X xY 
is compact) satisfies this condition. Indeed, if v is T-invariant, then: 

J f(x) dv(x,y) = J f(x,y) du(x,y) = 
= / f(T(x,y))du(x,y) = / f((T 1 (x,y),T 2 (x,y)))dis(x,y) = 



A similar reasoning can be applied to g. 

Given a continuous function c : X x Y — > [0,+oo), what can be said 
about 

I? 



a(c) := inf{ J cdrr : tt G n(T)}: 



What are the properties of optimal plans? We are interested here in 
Kantorovich Duality type of results. 
We will show the following: 

Theorem 4 (Kantorovich duality). a(c) is the supremum of the numbers 
a such that there exists continuous functions ip : X — > W, tp : Y — > M. 
satisfying: 

a + (p(x) - (piT^x, y)) + i/>(y) - i>(T 2 (x, y)) < c(x, y), W(x, y) e X x Y. 

We can list different interesting cases where we can apply the above result: 

1) If T doesn't depends of y G Y and T 2 doesn't depends of x G X, then 
we have the expression: 

a + <p(x) - tpfrix)) + il>(y) - 4>(T 2 (y)) < c(x,y) V(x,y) G X x Y. 
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In this case we are considering two variable invariant probabilities (one 
for Ti and the other for T 2 ) as marginals of an admissible plan. 

2) If X and Y are the Bernoulli space {l,2.,,<i} N , T\ = a is the shift 
acting on the variable x, (doesn't depend of y E Y) and T 2 = r x (y) (where 
Tj, j = 1, 2, d, are the inverse branches of a acting on the variable y) we 
have that T = a is the shift on {1, 2., , d} z and the above expression can be 
written as: 

a + (p(x) - (fi(cr(x)) + ip(y) - ip(r x (y)) < c(x, y) V(ar, y) E X x Y. 

In this case invariant probabilities 7r for the shift a : {1, 2., , <i} z — > {1,2.,, d} z 
are admissible plans. But not all admissible plan is a invariant. 

It is necessary to assume some special properties on c in order to get an 
optimal plan which is cr- invariant. This is the purpose of the next example. 

Example 5. Consider the points xq = (01)°°,xi = (10)°°, yo = (001)°°,?/i = 
(010)°°,?/2 = (100)°°. Let c a Lipschitz continuous cost satisfying: 

1) c(x ,y ) = c(xi,yx) = c(x ,y 2 ) = c(x 1 ,y 2 ) = 0, 

2) c > on the other points. 

It's easy to see that the unique optimal plan is given by: 

X r !r !r X r 

n = g<W?/o) + ^ d (xi,yi) + Q d (x ,V2) + Q d (xi,y 2 )- 

The support of this plan does not contain a a-invariant probability. 

Which is the right assumption to get an optimal plan tt which is a- 
invariant? In [8], [13], [H], [19], [20] some results in this direction are pre- 
sented considering a cost which is dynamically defined. 

Another class of examples: 

3) If X and Y are the Bernoulli space {1, 2., , d} N , and, for all (x, y) G X x 
Y, we have T x {x,y) = x, and T 2 (x,y) = r Xo (y), x = {x ,xi, ...) E {1,2., ,d} N , 
(where Tj, j = 1, 2, d, are the inverse branches of a), then, there is no (p(x) 
in this case, and, the above expression can be written as: 

a + ip(y) - ip{r x {y)) < c(x, y) V(z, y) E X x Y. 

This is the holonomic setting of [TT]. A duality result is proved in section 2 
in [TT] . Note that the ?/-marginal of a holonomic probability is invariant for 
the shift a acting on the variable y (see section 1 in [UJ). In the case c is 
Holder it is shown the existence of Holder realizers ip in sections 3 and 4 [TT]. 
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We will show in sections H] and E] here that the optimal plans can be ap- 
proximated by convex combination of optimal plans (of the classical transport 
problem) associated to measures supported in periodic orbits. In this way 
one can get an approximation scheme to the optimal plan based on a finite 
set of conditions. Our approach here is inspired in the point of view of tak- 
ing the temperature going to zero for Gibbs states at positive temperature 
in order to get results in Ergodic Optimization [21 J . The problem of fast 
approximation of maximizing probabilities by measures on periodic orbits 
plays an important role in Ergodic Optimization [TJ] [5] [7J. 

A paper which consider Ergodic Transport problems under a continuous 
time setting is [16J. 

We would like to thanks N. Gigli for very helpful comments and advice 
during the preparation of this manuscript . 

2 The case of one fixed probability and an- 
other variable invariant one 

Here we will prove Theorem [U We will adapt the main reasoning described 
in [27]. 

Given a normed Banach space E we denote by E* the dual space con- 
taining the bounded linear functionals from E to R. 
We will need the following classical result [27] [25] . 

Theorem 6 (Fenchel-Rockafellar duality). Suppose E is a vector normed 
space, E* its topological dual, O and 5 two convex functions defined on E 

taking values in R U {+00} . Denote 0* and S* ; respectively, the Legendre- 

Fenchel transform of O and S. Suppose there exists xq € E, such that 
xq G -D(O) D D{E), and, that O is continuous on xq. 
Then, 

inf [©(*) + S(x)] = sup [-e*(-/) - -*(/)] (4) 

x€E f eE * 

Moreover, the supremum in (HJ) is attained in at least one element in E*. 

Proof of TheoremUl First we prove I). 

We want to use Fenchel-Rockafellar duality in the proof. 
Define 

E = C(X xY)x MiY) 

where C (X x Y) is the set of all continuous functions inlxF taking values 
in R, with the usual sup norm || . W^. Moreover, M(Y) is the set of bounded 



7 



linear operators in C(Y) taking values in R with the total variation norm. 
Let P a (y ) be the set of invariant probabilities on Y. 
Define O : E — {+00} by 



9{u,v) = 



0, if u(x,y) > —c(x,y), V(x,y) G X x Y, and < 2 



+00, in the other case, 
Note that is convex. 

Define S : E — >RU {+00} by 



j'x^dn, if = <p(x) +ip(y) - (^oa)(y), 

with (tp,ij>) G C{X) x C(F),andz/ G P CT (F) 

+00, in the other case. 



Note that S is well defined. Indeed, if u — tpi(x) + V'i(y) — (V'i °")(y) = 
( fi2(x) +^2(2/) — (^2 then, integrating under any invariant probability 

f G P CT (y), we have that (pi(x) = (p2(x). Also note that S is convex. 
Observe that if 1/ G P a (Y), then (1, z/) G -D(O) H -D(S) and is continuous 
in 

Observe that 

inf [6(u, 1/) + S(u, 1/)] 
(u,u)eE 

= inf{ / ^ : + - o <r)](y) > -c(x,y), (y?,^) G C(X) x C(y)} 
Jx 

= m£{- f <pdw cp(x) + [i;-(i;oa)}(y)<c(x,y),( V >,i;)eC(X)xC(Y)} 
Jx 

= -sup{/ (pdn; <p{x) + [i>-{i>oa)}(y)<c(x,y),((p,4>)eC(X)xC{Y)} 
Jx 



= — sup / <f dji. 

{<p,i/j)e®c J 

Now we will compute the Legendre-Fenchel transform of and S, ini- 
tially, for any (n,g) EE*: by the definition of we get 

= SUp {< (-7T, I/) > -6(«,I/)} 

(u,v)eE 

= sup {-7r(u(a;,y)) -0(1/) : -u(ar, y) < c(x, y), \\v\\ < 2 } 
= sup {n(u(x,y)) - g(v) : u(x,y) < c(x,y), \\v\\ < 2}. 
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Following [27] we note that if ir is not a positive functional, then, there 
exists a function v < 0, v G C(X x F), such that, tt(v) > 0, therefore, taking 
u = Xv (remember that c > 0), and considering A — > +oo, we get that 

sup {ti(u) : u(x, y) < c(x, y)} = +oo. 

When 7r G M + (X x K) and c G C(X x K) we have that the supremum of 
Ti (u) is, evidently, 7r(c). 
Therefore, 

tt(c) + sup -^(i/), if 7r G M + (X x y) 

e*((-7r, -«/)) = < IHI - 2 (5) 

+oo, in the other case. 

Analogously, by the definition of H we get that 

S*(tt,^)= sup {< (7r,g),(u,u) > -E(u,u)} 
(u,v)eE 

{n(u(x,y)) - f ipd/JL + g(u) : 
y) = (p(x) + ip(y) — ip(a(y)) where (<p, ip) E C{x) x C(Y) 



sup <.7c((p(x)+t/)(y)-i/)(<T(y)))- <pdn + g(v) 

(<p,iP)eC(x)xO(Y), veP a (Y) I ./ 

If 7r(y(a;)) ^ J Lpdfi (we can suppose greater) for some (p, taking X.ip 
and A — > oo, the supremum will be equal to +oo. Analogously if 7r(ip(y) — 
ip(cr(y))) 7^ (we can suppose greater) for some ip, taking X.ip and A — > oo, 
the supremum will be +oo. 

In order to simplify the notation, define 

ir((p(x)) = fipdn and 7r(ip(y) - ip(a(y))) = 

v(tp, vo g c(x) x c(y) 



rr(/i) = |vr g m(a x y) 

We just show that 



7T, 



sup g(u), if 7r G 

(6) 

+oo, in the other case. 

We know that the left hand side (jlj) is equal to — sup^ ^) £ $ c / y? ^A*, and also 
by © and we know that the right hand side of (jjj) coincide with 
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sup 

(T,ff)£E* 



sup 

(77,g)eE* 



tt(c) + sup -g(u) + sup g(v), if tt G M + (X x Y) n II*(/x) \ 

||i/||<2 ^eP CT (y) 

+00 in the other case 



-tt(c) + inf g(v) - sup c/(z/), if 7r G M + (X x Y) D IT(» 
IMI<2 ^eP CT (y) 



in the other case 



= sup{-7r(c), tt G M+(X x Y) n II*(/i)}, 
where the last equality is obtained taking g = because ||z/|| = 1 for any 

^ e P CT (r). 

Finally, we observe that if 7r G M + (X x F) fi II* (/i, z/) then: 

7T(1) = /i(l) = 1, 

7r(u) > 0, if M > 0, 
7r is linear. 

From these properties we get that 71 G P(X x K). 

Moreover, by definition of IT(/i), the projection of tt in the first coordinate 
is fi, and, the projection of it is invariant in the second coordinate. It follows 
that M+(X xY)n n*(/x, 1/) = n(/x, a). 

Therefore, from this together with (j3J) we get 



sup / if d\i = — inf / c dir 



sup (fd/i— inf / cc?7r. 
((p,il>)e& c J 7ren(M,o-) 7 

Note that theorem [6] claims that 

sup [-9*(-/)-H* (/)] = - inf / cdrr 

is attained, for at least one element, and this shows the existence of an 
optimal plan. This shows I). 

After we get the probability v we can consider the classical transport 
problem for fi, v and c, and, finally, we can get some well known proper- 
ties described in the classical literature (as, slackness condition, c-cyclical 
monotonicity, etc.). 
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Now, we will prove II). This will follow from the following claim. 

Claim: Let X be a compact metric space, Y = {1, d} N , c:IxF-^l 
be a Lipschitz continuous function and /jbea probability measure acting in 
X. Let 7r G II(/i, a) minimizing the integral of c. Then, there exist Lipschitz 
continuous functions ip(x),ip(y) such that: 

i) <p(x) + ip(a(y)) - $(y) < c(x, y) 

ii) f (fi(x) d\i = f cdn. 

Let (3 the Lipschitz constant of c. 
First note that given continuous functions ip and ip satisfying 

<p(x) + ip((r(y)) - ip(y) < c(x, y), 

then, there exists Tp and ip, Lipschitz functions with Lipschitz constant j3 
satisfying: 

Tp(x) + ip(<r(y)) - ip(y) < c(x, y), 

and, 

Tp > V- 

We can choose tp satisfying < ip < (5. 

Indeed, for any a n (w) = y and Xq, ...,x n _i G X: 

n-1 

This shows that 

n-1 

ip(y) := inf{^ c(x h a l {w)) - tpfa) : n > 0, a n {w) = y, x { G X} 
is well defined. 

We remark that ip is a Lipschitz function with the same constant (5. Note 
also that 

</?(x) + ip(a(y)) - ijj(y) < c(x, y). 
Now for each x fixed, define Tp{x) as the greatest number such that for any 

■ ' + ip(v(y)) - ip{y) < c(x, y). 

We note that Tp > tp and Tp(x) = inf{c(:r, y) + ip(y) — ip((r(y))}. We also note 

y 

that Tp is a Lipschitz function with same constant (3. We remark that we can 
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add a constant to ip, and, so we can suppose without lost of generality that 
< i) < (3. 

Now, we prove the main claim. 

By (jSJ), there exists sequences of continuous functions (p n and if) n , n e N, 
such that 

ifn(x) + 1p n (v(y)) - 4>n(y) < c(x, y) 

and 

(p n djJ, — > I CdlT. 



From the above reasoning we can get <p n ,i[) n Lipschitz continuous func- 
tions such that 

9 n { x ) + - $ n (y) < c ( x > y)> 

and, 

Tp n djj, — > I cdir. 



For a fixed e > 0, we get 

cdn — e < I Tp n dfx < I c dir, 



for n large enough. Particularly, for a fixed n sufficiently large, there exist 
x n , x' n G X such that 

cdn — e < Tp n (x n ) and Tp n (x' n ) < J cdir. 

Using the fact that Tp n is Lipschitz continuous, with constant (3, and, denoting 
by D the diameter of X, we conclude that for large n: 

cdn - e- D (3 <Tp n < Jcdn + Df3. 

So we can apply the Arzela-Ascoli theorem and, finally, we get continuous 
functions <p, ip satisfying: 

i) tp(x) + ip(a(y)) - ip(y) < c(x, y) 

ii) f f(x)dfi = f cdn. 

We know from the first reasoning that we can assume ip and ip are Lips- 
chitz continuous functions. This shows II). 

□ 
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3 Generic properties: a unique optimal plan 



Lemma 7. Let K be a compact set in M 2 , and for each r > 0, define K r as 
the set of points (x, y) G K such that x + ry is maximal. Then de diameter 
of K r converges to zero when r — > 0. 

Proof. See [3] page 306 for the proof. □ 

Corollary 8. With the hypothesis of the above lemma, for each e > 0, there 
exists a r > 0, such that, for r > r > and (xx,yi), (x2, y-i) G K r , we have: 

1 2/i - 2/2 1 < e. 

The bellow theorem follows from the same arguments used in proposition 
9 in [3]. 

Theorem 9. Let X be a compact metric space, Y = {l,...,d} and \x be 
a probability measure in X. Let C(X,Y) be the set of continuous functions 
from X xY to IR + with the uniform norm. The set of functions c G C(X, Y) 
with a unique Optimal Plan in IT(/i, a) is generic in C(X,Y). The same is 
true for the Banach space H(X, Y) of the Lipschitz functions with the usual 
norm. 

Proof. On this proof, we are going to consider tt an optimal plan if J cdn is 
maximal (just consider the change of c by — c). 

We start studying the space C(X,Y). Given an countable family (ej) igN 
dense in C(X, Y), the set of functions in C(X, Y) with two or more optimal 
plans coincides with: 

m,n£N 

where 

X m ,n '■= { c £ C(X,Y) : 3 7r, x optimal plans, / e n d(n — x) > — }• 



m 

Then it's sufficiently to prove that X m ^ n is a closed set with empty interior. 

Claim 1: X m , n is a closed set. 
Indeed, we note that C(X,Y) is a normed space. Consider c s in X m ^ n con- 
verging to c (when s — > oo). Let (tt s , x s ) be the optimal ones associated with 
c s in X m n . We can suppose, by taking a subsequence, that n s — > ir, Xs — > X-, 
where tt, x are probability measures in X x Y . So 

e n d(n - x) = lim / e n d(ix s - x s ) > — ■ 



s->oo / 777, 
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Clearly vr, x G H(fi,a). Also, by the above relation, they are different mea- 
sures. We want to show that the limit function c is in X m n . We only need to 
prove that it and x are optimal plans to c. Suppose by contradiction there 
exists ( G a) such that f cd( > f cdn + e. So for s large we have: 

c s dC,> Jcd(^ — e/3> Jcdn + 2e/3> J cdn s + e/3> J c s dir s . 

This is impossible because n s is an optimal plan for c s . Therefore, \ is an 
optimal plan for c. 

Claim 2: X m n has empty interior. 
Indeed, for a fixed c G X m ^ n we can show that c + re n ^ X m>n when r > is 
sufficiently small. Consider 

K = | (^J cdn, J e n d-n^j : n G Il(/i, a 

if is compact and contained in M 2 . Then, by the Corollary [HJ when e = l/2m, 
there exist a r , such that for r > r > we get: if f(c + re n ) dn and 
/(c + re„) g?x are maximal (this means tt, x optimal plans to / + re n ), then 
|J (e n ) d(n — x)\ < ^ = l/2m. This show that c + re n ^ X mi „. 

In the space H(X, Y) we can get similar results. This can be obtained 
with the same arguments used before together with the following remarks: 

a) A dense enumerable family {e n } in H(X, Y) will be a dense sequence in 
C(X, Y). In this way two elements n,x G II(/z, a) will be different, if and 
only if, J e n d(n — \) 7^ 0, for some e n . 

b) Moreover, the set 

X m n := {c G -ff (X, 30 : 3 7r, x optimal plans, / e n d(ir — x) > — } 

J m 

will be a closed set by the same arguments used above. We can also show, in 
a similar way, that it has empty interior, but we note that we need c + re n G 
H(X, Y), and c + re n — > c in Lipschitz norm. This is true because we can 
consider e n G H(X,Y). 

□ 



4 Zeta-measures and Transport 

When fi, v have support in a finite number of points, then, the optimal plan 
7r(^, v) for a cost c in the Classical Transport Theory can be explicitly ob- 
tained by Linear Algebra arguments [27]. Indeed, suppose 

/i = a x 8 Xl + ... + a n 5 XnJ 
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and, 

u = b 1 5 yi + ... + b m 5 ym . 

Then any transport plan n have support contained in {x±, ...x n }x{yi, ...,y m }. 
Denoting by 7Tjj the mass of n in (xi,yj) we have that the variables ir^ need 
satisfies the linear equations: 
vertical equations: 

0-1 = TTll + ••• + TTml 
CL n Tfln ~\~ ■■■ ~\~ T^mn 

horizontal equations: 

b\ = 7Tn + ... + 7T ln 
bm TTml ■ ■ ■ ^mri 

The set of solutions of this equations defines a convex in M ra - m . The 
conditions > will restrict the solutions to a bounded convex set with 
finite vertexes. So given a cost function c : X x Y — > [0, oo), denoting their 
restricted values by Qj := c(xi,yj), we have that the optimal plans to (/i, v) 
are the point of the above convex set that minimize the linear functional: 

^ \ C ij^ij- 

By convexity arguments there is an optimal point in the vertexes of the 
underlying convex set. The conclusion is that by Linear Algebra arguments 
we can find a finite number of points such that at least one of these will be 
optimal for the integral of c with the given marginals /i and v. 

Note that these finite vertexes points are determined before we consider 
the given cost function. 

It is natural in the Ergodic Theory setting to try to approximate a general 
invariant probability by the ones which posses the simplest behavior: the 
periodic probabilities. These are the ones that we can make computations 
more easily. 

We note that to minimize the integral of the cost function c is the same 
that to maximize the integral of the function — c. The plan that realizes 
this optimal integral will the same if we add a constant to — c. Bellow we 
consider the problem of finding a transport plan maximizing the integral of 
a cost c strictly greater than zero. A transport plan from // to v maximizing 
the integral of c will be called a maximizing plan. Below we consider a 
compact metric space X and Y = {!,..., d} N . 
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Definition 10. For fixed \x G P(X) and a continuous function c : XxY — >• M. 
we define a probability measure in X x Y by the linear functional n : 
C(X x 7) 4 1, suc/i i/iai, /or eac/i u> G C(X x Y) associate the number: 

J2ueFi Xn eP ' n ' Ic{x ' v)d7rM I wdir(ji,v) 

Epfi.n. J c(x,y)dn(jj,,v) ' 
p£Fix n ° 

where Fix n denotes the set of invariant measures in Y supported in a periodic 
orbit of length n, and, 7r(/x, v) denotes a maximizing plan from /i to v with 
cost function c (we don't impose other conditions on the chosen the plan). 

In the case \i is supported in a unique point xo, we can define the function 
A(y) = c(xq, y), and this measure can be written as: 

V c 0.A"(y) w n (y) 
/^y£Fix n ° 

where w(y) = w(x ,y) and A n (z) = A(z) + ... + A(a n ~ 1 (z)). This kind of 
measure (also called zeta measure) is considered in Thermodynamical For- 
malism [24] and they can be used to approximate Gibbs states, and, also 
the measure that maximizes the integral of A among the invariant measures 
(see |21j). Therefore, in some sense, the above defined family of probabilities 
extend a well known concept used in Ergodic Optimization. 

In the case of that fi have finite support these zeta-measures can be 
determined by Linear Algebra arguments like we remarked above. 

Remember that II(/i, a) is the set of probabilities measures that coincides 
with fi in the first coordinate and is invariant in the second coordinate. 
The next theorem follows the ideas used in thermodynamical limit when 
P,n-> oo |2I]. 

Theorem 11. When (3, n goes to infinite, any limit measure tt^ of convergent 
subsequence of (p^, in the weak* topology, belongs to H(fj,,a). Moreover, if 
c > 0, then, ir^ maximizes the integral of c among the measures in U(fi,a). 

Proof. We begin by proving that for 0, n fixed, the corresponding zeta- 
measure is in Tl(fi, o~). Let w a function depending only on x. Then 



Epfi.n. f c(x,y)dn(fj,,u) 
u£Fix n c 

_ Y.„ &Fl xy- n - Sc{x ^ M 

Epfi.n. f c(x,y)dn(fj,,i/) 
v£Fix n c 



16 



w d/j, 



u&Fix n ° 



E p j3.n. J c(x,y)dn(fi,u) 
ueFix n ° 

W dfX. 

Now, consider a fixed function w depending only of y. Then, we have: 

(,p,n(WO(T)— J c(x,y)<hr(u,v) 

Z^iueFix n c 

y e P-n. J c{x,y)dK(n,v) f WO(TC lu 

tr T>Xn u 

Ep/3.n. f c(x,y)dir(fj,,w) 
ueFix n c 

T,yeFix n eP ' n ' !<X,y)M>l ' V) I wdu 

Epfi.n. J c(x,y)dn(fj,,u) 
u£Fix n c 

E„gF*r„ e ^ n ' 1 c{x ' y)dnM J w dn(n, V) _ 

V- e f3.n.Jc(x,y)dK{n,v) <»PA W )- 

Z^v&Fix n ° 

This show that £p >n G n(/x, cr). So, when /3,n goes to infinite, any limit 
measure of convergent subsequence of C/3,n is on IT(/i, a). 

Suppose C^^j ->■ TToo, when j -)> oo. 

Let 7r* G a) maximizing the integral of c. Let i/* the projection 
of 7r* on the second coordinate y. Then, v* is a invariant measure. Let 
ZV £ Fix n . be a subsequence converging to i/* in the weak* topology. If 
7r ra . is a maximizing plan from to u n ., then, there exist a subsequence 7r nj 
converging to a maximizing plan ir from /i to v* ([28J page 77). It is easy to 
see that 

cdn = I cdrr* 



and, therefore, 7r is maximal. In other words, it maximizes the integral of c 
among the measures in IT(/i, a). We denote this integral by 1(c). We want 
to show that ^(c) > /(c), where the subsequence CPi,m converges to in 
the weak* topology. From the above arguments we know that: 



Given e > 0, for sufficiently large i there exist v G Fix ni such that: 

J cdn(fi, v) > 1(c) — e. 
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Take e > 0, such that (/(c) — e) > 0, and define: 

A nX e ) = i v e Fix ni ■ J cdn(ii, v) < 1(c) - e} 

BnX e ) = {v e Fix n% : J cdn(fj,, v) > 1(c) - e}. 



Then, we have: 

e Pi-ni. J c{x,y)d-K{n,v) < e ft.ni.(/(c)-e) 
veA ni {e) veA ni {e) 

< gnilogCdHft-ni-UW-e) 

and 

£ e P i - ni -S'*- x ' v) ' hc ^ v) ( cdix(^v) < e niXog{d)+ ^ n ^ I{c) ~ £) (I(c) - e). 

By other hand, if is sufficiently large, B ni (e/2) is non empty It follows 
that 

g/Si-fii. / c(x,y)dn(fi,v) y ^.n,. J c(x,y)dn(^u) 

f&B ni {e) ueB ni (e/2) 
> e f3 t .n t .(I(c)-s/2) ^ 

and, 

£ ePi-m-S«*>v)**M J c dn(}i,v) > e^ n '- {I(c) - £/2) (I(c) - e/2). 



Then, we get 



„ „ ,■ ^GA ni ( £ ) e 

U < hm — 5 — — — — - 

woo V „ , , e Pi-iH- J c(x,y)dir(ix,v) 



e m log(d)+/3j.7ij. (J(c)-e) 



Moreover, 



< lim .„ . = lim e ^^)-p t .n t .e/2 = Q 

i— >oo QPi-ni.(l(c)—£/z) j— >oo 



, E^„.( £) e ft - ri! - /c( ^ ( ^ ) /crf7r( /U ,z/) 

< hm ! < 

" i4M E, e j)„.( £ ) eft '"" /cMilM fcdn(fx,u) ~ 

e n i log(d)+/3 i .n i .(I(c)-e)^j^ _ £ ^ 



e /W(/(c)- £ /2)(/( c ) _ e /2) 
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Finally, 



lim inf 



lim c "» lo K^)-ft-"»-£/ 2 J ( c ) e = o. 
i->-co /(c) — e/2 



= lim inf 5= 7, t-, — m — ; 

j-s-oo V „ . . e Pi-n%. ] e{x,y)dnt(jx,v) 

Z^ Ui eB n (e) c 

> lim inf { { ) } 

> /(c) - e. 

Taking e — > 0, we get: 

h S£ f y e ^. Ic (x, y)d , M ^ /(c) - 

Using the fact that Cft,n; — >• ^<x> we conclude that ^(c) > 1(c). □ 

5 Two variable invariant probabilities and other 

cases 

We start this section by proving Theorem HI The proof follows basically the 
same kind of ideas that were used in Theorem [TJ 

Proof. Define 

E = C(X xY)x M(X x Y) 

where M(X x Y) is the set of bounded linear functionals from C(X x Y) to 
R with the norm given by the total variation. 
Let : E — > R U {+00} given by 



0, if u(x, y) > —c(x, y), V(x, y) G X x Y, and \\u\\ < 2 

+00, in the other case, 
Note that G is a convex function. 
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Define E : E — > R U {+00} by 



a, if u(x,y) = a + tp(x) - v{Tx{x,y)) + ijj{y) - if>(T 2 {x,y)), 
with (<£>,V) e x C(r),andz/ G n(T) 

+00, in the other case. 



We remark that H is a well defined convex function. 

If ^ G n(T), then (1, v) E D(Q) fl and 6 is continuous in (1, v). 



Note that: 



inf [6(w, z/) + H(it, v)\ 



inf r ](^) 6 c(i)xc(n 1 

\ ' a + y?(x) - ip{Ti{x,y)) +ip(y) -if>(T 2 (x,y)) > -c(x,y) J 
^i)GC(I)xC(F), 



ml a ' a + <p(z) - <f(Tx(x,y)) + ^(y) - ip(T 2 {x,y)) < c(x,y) 

sup<«- G C(X) x C(y), 

a + <p(z) - <p(Ti(x, y)) + ^(y) - ^(T 2 (x, y)) < c(x, y) 

So the left size of (jlj) is: 

-sup la- 6 C7(X) x C(y), 1 

[ a + </?(x) - y{T x (x,y)) + i/;{y) - i/;(T 2 (x,y)) < c(x,y) J 

Now we will compute the Legendre-Fenchel transform of and 5. For 
each (-7T, g) E E*\ 

©*((- 7r 7 ~g)) = sup {< (-n,-g),(u, v) > -6(u, 1/)} 

= sup {-7r(u(x,y)) -(/(z/) : y) < c(x, y), ||z/|| < 2 } 

(u,v)eE 

= sup M«(a:,y))-^(i/) : u(ar, y) < c(x, y), \\u\\ < 2} . 

Following [27] note that if tt ^ M + (X x K), then there exists a function 
t> < in C(X x K), such that, 7r(v) > 0, so taking u = Xv (remember that 
c > 0), when A — >■ +00, we have that sup ueE {tt(u) : u(x,y) < c(x,y)} = 
+00. 

Moreover, if ir E M + (X x Y), using the fact that c E C(X xY), then we 
have that the maximum of tt(u) is given by 7r(c). 
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Therefore, 

e*((-7r,-s)) = 

Now we analyze 5*: 



tt(c) + sup -(/(i/), if vr G M+(X x F) 
IMI< 2 

+00, in the other case. 



E*(n,g)= sup {< (vr, g), (u, v) > -E(u, u)} 

{u,L>)eE 

ir(u(x,y)) -a + g(u) : 
sup ^ u(x,y) = a + <p(x) - ip(Ti(x,y)) + ip(y) - i/j(T 2 (x,y)), 
(«.")eB [ W ith 0,-0) e C(a;) x C(Y")andz/ G n(T) 

sup 7t(q;) — a + 7r((y9 — o T x ) + 71- (-0 — ^ o T 2 ) + g(z^). 

a,((p,v)6C(x)xC(y),i/en(T) 

If ir(<p(x)—<p(Ti(x, y))) 7^ (we can suppose greater than zero) for some <p, 
by taking X.ip and A — > 00, the supremum will be equal to +00. Analogously, 
if Tr(il)(y) — tp(T>i(x, y))) 7^ 0, the supremum will be +00. If 7r(l) 7^ 1 (we can 
suppose greater than one), then taking a — > 00, the supremum will be +00. 

Define 

7r(l) = 1, 7r(v? - o T x ) = 7r(^ - ^ o T 2 ) = 0, 
V(^,f)6C(X)xC(F) 



n*(T) = s 7r G M(X x F) : 



Therefore, 



sup{^) : v G n(T)}, ifvrGn*(T), 
+00, in the other case. 



(9) 



We know that the left size of PJ is given by dU) . By ([S]) and © , the right 
size of (HI) will be: 



sup 

(7T, 9 )eE* 



/ (tt(c) + sup{-g(v) : < 2}) + sup{g(v) : 1/ G II(T)}, \ 
if vr G M + (X x Y) fl n*(T) 



y +00, in the other case 



/ 



sup 

(7T, 9 )eE* 



/ (-tt(c) + mi{g(u) : H < 2}) - sup{g(v) : 1/ G U(T)}, \ 
if 7r G M+(X x y) nn*(T) 



y —00, in the other case 
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/ 



= sup{-7r(c), 7T G M + (X xY)n IT(T)}, 
where the last equality is obtained taking g = 0. 

We remark that if vr G M+(X x Y) n II* (T) then: 
?r(l) = 1, (by definition of IT(T)) 
7r(u) > if u > 0, 
7r is linear in C(X x K), 

Therefore, we have 7r G P{X x Y), and, by definition of II* (T), we get that 
vr G IT(T). The conclusion is that M + (X x Y) nn*(T) = n(T). So the right 
side of (j4]) will be: 

sup{-vr(c), vr G II(T)} = -inf{7r(c), vr G II(T)}. 

Therefore, we conclude from (HI) that: 



3(^)6C(I)xC(n 

a + tp(x) -ip(Tx(x,y)) + ip(y) - ip{T 2 (x,y)) < c(x,y) 

= -inf{7T( C ), 7T G II(T)}, 

or, in another form that 

inf{7r(c), 7T G II(T)}. 



' " !) ' a + y?(£) - <p(Ti(x,y)) + i>(y) -ip{T 2 (x,y)) < c(x,y) 



□ 

Proposition 12. Suppose c is continuous. Denote a = inf / c(x,y)dir. 
// i/iere ezisi y> G C(X) and ^ G C(Y) satisfying 

a + ip{x) - ip(T x {x)) + 4>{y) - 4>{T 2 {y)) < c(x, y) V(x, y) G X x Y, (10) 

inf {c(x, y) + ... + c(T n_1 (x, y)) — net : n > 1 and(x, y) E X x Y} > — oo. 
Proof. Suppose by contradiction that 

inf{c(s, £/) + ...+ c(T n_1 (x, y)) — na : n > 1 and (x, y) G X x Y} = — oo. 

Also suppose that there exists <p and ip satisfying fflUj) . Then we have 

inf{^(x)-^(T 1 "(a;,y))+V(y)-V(T 2 ri (x, 2 /)) : n > 0and(x,y) G XxY} = -oo. 

This is impossible because X and Y are compact sets and (p and ip are 
continuous functions. □ 
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Proposition 13. Suppose X = Y = {0, 1} N ; T X =T 2 = a, and, that c is a 
Lipschitz function, then there exists (p(x) andip(y) Lipschitz continuous such 
that 

a + <f(cr(x)) - tp(x) + ip(cr(y)) - tp(y) < c(x, y). 

Proof. Denote by (3 a Lipschitz constant for c. 

By definition of a there exist a increasing sequence a n — > a and continu- 
ous functions <p n ,ip n such that: 

a„ + (p n (a(x)) - ip n (x) + ip n (cr(y)) - ipn{y) < c{x, y). 

From this relation we have that if a rn (z) = x and y , ...,y m _i belongs to Y, 
then 

m— 1 
i=0 

Therefore, 

m— 1 

M{^2c(a l z,yi)+ip n (yi)-ip n (a(yi))-a n : m > 0, cr m (z) = x, yi E Y} > -oo. 

i=0 

Denote by ^ n (x) this infimum. Note that this function satisfies: 
a n + Tp n {a{x)) - 7p n (x) + i/> n (<r(y)) - ip n (y) < c(x, y). 

It is easy to see that Tp n is Lipschitz continuous with the same Lipschitz 
constant /3 of c. Using this last inequality and the same arguments used 
before (applied to ip) we can construct a Lipschitz continuous function ip n 
with the same Lipschitz constant (3 satisfying: 

"« + V n {?(x)) - lp n (x) + tp n (a(y)) - i) n {y) < c{x, y). 

Note that we can add constants on Tp n and ip n and the conclusions are the 
same. Then, we get: there exists Lipschitz continuous functions Tp n and ijj n 
with Lipschitz constant j3, which bounded bellow and above, respectively, by 
and j3, such that: 

Now using Arzela-Ascoli theorem we obtain continuous functions ip and ip 
satisfying: 

a + <p(<r(x)) - tp(x) + ip(<r(y)) - ip(y) < c{x, y). 

Applying the same reasoning of the previous arguments we can also construct 
Lipschitz functions satisfying this inequality. □ 
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Proposition 14. Let C(X, Y) be the set of continuous functions from 1x7 
to M + with the uniform norm. The set of functions c £ C(X,Y) with a 
unique Optimal Plan in Tl(T) is generic in C(X,Y). The same is true for 
the Banach space H(X,Y) of the Lipschitz functions with the usual norm. 

Proof. The result follows from adapting the proof of Theorem |9j □ 



6 Zeta-measures for the second class of prob- 
lems 

In this section we suppose X = Y = {0, 1} N and T% = T 2 = a is the shift. 
On this case we have II(T) = II(cr) is the set of probabilities 7r in X x Y such 
that project on er-invariant measures in X and Y. 

Bellow we consider the problem of finding a transport plan in U.(a) max- 
imizing the integral of a cost c strictly greater than zero. A transport plan 
maximizing this integral will be called a maximizing plan. By changing 
the signal of the cost we can get from this the analysis of usual minimization 
problem of Transport Theory. 

Definition 15. For a fixed cost c we define a probability measure in X xY by 
the linear functional Qp yn : C(X x Y) — > R ; such that, to each w £ C(X x Y) 
we associate the number: 

Epfi.n.j c(x,y)dn(fj,,u) 
fj,,uGFix n 

where Fix n denote the set of invariant measures in X = Y supported in a 
periodic orbit of length n, and, 7r(/x, v) denote a maximizing plan from \l to 
v with cost function c (we don't impose other conditions on the plan). 

This zeta-measures can be determined by Linear Algebra arguments. In- 
deed, note that if v £ Fix n , then the plan 7r(/x, v) can be determined by 
the study of certain permutations (see page 5 in [27J). 

Theorem 16. When (3, n goes to infinite, any limit measure -k^ of convergent 
subsequence o/C/3,«; ^ n the weak* topology, is on n(cr). Also, if c > 0, then, 
iToo maximizes the integral of c between the measures in n(cr). 

Proof. We begin by proving that for /3, n fixed, the zeta-measure is in II(<t). 
Let w be a function depending only on y. Then 



[w o a) 



Epfi.n. J c(x,y)dn(tj,,u) 
li,u£Fix„ 
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E e f3.n. f c(x,y)dir(fi,u) 
p p.n. $ c(x,y)dir(n,v) f ... J y 

Epfi.n. f c(x,y)dir(fi,u) 

^-^fi,u£Fix n 

If w depends only on x the argument is similar. This shows that C/3,n G n(cr). 
Then, when /3, n goes to infinite, via a convergent subsequence, any limit 
measure it^ of (p jTl , in the weak* topology, will be on n(er). 

Suppose C/3j,nj ->■ TToo, when j -)> oo. 

Let 7r* G n(cr) be a probability maximizing the integral of c. Let 
respectively, the projection of tt* in the first and second coordinates. Then, 
fjL*,u* are invariant measures. Let /x n . , i/ n . G -Fix nj . subsequences converging 
to /x*, v* in the weak* topology. If 7r nj is a maximizing plan from \i n . to z/„ 3 ., 
then, there exist a subsequence 7r ni converging to a maximizing plan tt from 
//* to v* ([25] page 77). Clearly 



Cd7l = / CG?7T 



and, therefore, 7r is maximal. This means that ir maximizes the integral of 
c among the measures in n(er). We denote this integral by /(c). We want 
to show that r K o{c) > I( c )- We note that the subsequence Cft,™, converges to 
7TQO in the weak* topology. From the above arguments we get: 

given e > 0, for sufficiently large i there exists u., v G Fix nv such that, 

J cdir(u,, v) > 1(c) — e. 

Consider a fixed e > 0, such that, (1(c) — e) > 0, and, define: 

4*i(e) = {(a*, f) e Fzx nj : / cck(/i, »/) < 1(c) - e}, 



B ni ( £ ) = ^) e : J cdn(/j,, v) > 1(c) - e}. 

Then, we have that 

gft-"i- / c(x,y)dw(n,u) ^ e (5i.ni.(I(c)-e) 
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< e 2n i log(2)+/3 i .n i .(/(c)-e) 

and, 

J2 e ^ n ^ c{x ' y)dnM I Cdn(n, V) < e 2n i log(2)+/3 i .n i .(7(c)- £ )^^^ _ £ y 

By other hand, if is sufficiently large, B ni {e/2) is not empty. Moreover, 

gft-"i- / c{x,y)d%{n,v) y gft-ni. / c{x,y)d-K{n,v) 

{li,v)eB ni (e) (p,u)eB ni (e/2) 

> e Pi. ni .{I{c)-e/2) 

and, 

£ eft^-M**)*^.") / crfvr(/i, i/) > e ft-««-W C )-^/2)(j( c ) _ e /2). 

Then, 

0<lim^ )6A - (£) 



ppi.m. f c(x,y)dir(ji,i>) 



and, 



2n; log(2)+/3 i .n i .(7(c)-£) 

< lim 5 ,„ . = lim e 2™ I log(2)-ft.n l . £ /2 = q 

~~ i^oo e Pi.ni.(I(c)-e/2) 



U < hni — - t—. — — r— - 

~ *■+<*> J2 M eB ni(£ ) e^-M^M Jcdn(fx, v) 

e 2n i log(2)+ft.n i .(/(c)-£)^j^ c ^ _ £ ) 
~ /™ e ft.n i .(7(c)- e /2)( / ( c ) _ J/2) 

= lim c 2 "» lo s( 2 )-ft-"»-£/ 2 ^ ~ £ = o. 
j-^oo /(c) — e/2 



Therefore, 



lim mt — — 1 t- — ——, — : 

i~>oo V, n „ pPi-ni. J c{x,y)dir(n,v) 



lim inf 



i-»oo V/ \ „ , . e Pi-ni. f c(x,y)dir(p,u) 

Z^i(n,u)eB n , (e) 



/(/i,i/)6B„ j (e) ' 
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* h £l nf y f \ R ( ,^.j c{x , y)d , M 

l~,{w)eB nt (e) c 

> 1(c) - e. 

Taking e — > we get 

liminf E^) effa!B< eft-"*-J^)^)/ C d7r(/x,i/) > 

i-»-oo V, , „ pPi.rii. j c(x,y)dTr{n,v) ~ ^ >' 

Then, using the fact that C/3 „ converges to n^, we finally get ^(c) > 
/(c). □ 
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